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ABSTRACT 

Context. The advent of space-borne photometers such as CoRoT and Kepler has opened up new fields in asteroseismology. This is 
especially true for red giants as only a few of these stars were known to oscillate with small amplitude, solar-like oscillations before 
the launch of CoRoT. 

Aims. The G6 giant HR2582 (HD 50890) was observed by CoRoT for approximately 55 days. We present here the analysis of its 
light curve and the characterisation of the star using different observables, such as its location in the Hertzsprung-Russell diagram and 
seismic observables. 

Methods. Mode frequencies are extracted from the observed Fourier spectrum of the light curve. Numerical stellar models are then 
computed to determine the characteristics of the star (mass, age, etc..) from the comparison with observational constraints. 
Results. We provide evidence for the presence of solar-like oscillations at low frequency, between 10 and 20yuHz, with a regular 
spacing of (1.7 ± 0.1)yuHz between consecutive radial orders. Only radial modes are clearly visible. From the models compatible with 
the observational constraints used here, We find that HR 2582 (HD 50890) is a massive star with a mass in the range (3- 5 M Q ), clearly 
above the red clump. It oscillates with rather low radial order (n = 5-12) modes. Its evolutionary stage cannot be determined with 
precision: the star could be on the ascending red giant branch (hydrogen shell burning) with an age of approximately 155 Myr or in 
a later phase (helium burning). In order to obtain a reasonable helium amount, the metallicity of the star must be quite subsolar. Our 
best models are obtained with a mixing length significantly smaller than that obtained for the Sun with the same physical description 
(except overshoot). The amount of core overshoot during the main-sequence phase is found to be mild, of the order of 0.1 H p . 
Conclusions. HR2582 (HD 50890) is an interesting case as only a few massive stars can be observed due to their rapid evolution 
compared to less massive red giants. HR 2582 (HD 50890) is also one of the few cases that can be used to validate the scaling relations 
for massive red giants stars and its sensitivity to the physics of the star. 

Key words. Stars: oscillations - Stars: individual: HR2582 (HD 50890), HD 50890, HIP 33243 



1. Introduction 

Solar-like oscillations (p modes excited by turbulent motion in 
the convective outer layers of the star) have first been observed in 
a couple of red-giant stars using either spectroscopic data from 
the ground (see e.g. Bedding & Kjeldsen 2007) or space pho- 
tometric data (Barb an et al.|2007| l. The detection of such oscil- 
lations in evolved stars, such as red giants, is made difficult by 
their very low frequencies (generally below 100 //Hz), but made 
easier by their intrinsically larger amplitudes (around 100 ppm 
in photometry) and smaller linewidths compared to solar-like 
stars. The use of long, uninterrupted observations provided by 
the space-based photometer CoRoT (Baglin et al. 20061) has per- 



mitted the detection of oscillations in hundreds of red giants (De 
Ridder et al. 2009). The observed spectra show different char- 



* The CoRoT space mission, launched on 2006 December 27, was 
developed and is operated by the CNES with participation of the 
Science Programs of ESA; ESA's RSSD, Austria, Belgium, Brazil, 
Germany and Spain. 



acteristics: regularly spaced frequencies or more complex spec- 
tra (|Hekker et al.|20 09). Analysing Kepler data, Bedding et al. 
( |2010| l found many of these regularly spaced frequency spectra. 
Based on the same mission results, the fundamental seismic pa- 
rameters of red giants were described by Kallinger et al. ( 2010a| l. 
Observed complex or regularly spaced spectra can be compared 
to the theoretical expectations of Dupret et al. ( 2009 1 and can 
be understood as the presence or the absence of mixed modes. 
The regular spectra can be described using a "universal pattern" 
(Mos ser et al.||2011b|l. The observed pulsating stars described 
in |Hekker et al.| (|2"009) show a distribution of seismic parame- 
ters that indicate that they belong to the so-called "red clump" 
of low mass, He-burning, evolved stars (Mi glio et al.|2 009). The 
determination of the mass, radius and age of these evolved stars 
is a challenge (Kallinger et al.|2010b) > when their luminosity and 
effective temperature are poorly estimated. Scaling relations de- 
scribing seismic characteristics are then used to determine their 
mass and radius. Few individual red-giant stars have been mod- 
elled using seismic constraints. For example, Miglio et al. ( 2010) 
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Fig. 1. Time series of the relative intensity variations of HR 2582 
(HD 50890) obtained by CoRoT. 



found the signature of sharp structure variation in the red gi- 
ant HR7349 observed by CoRoT and analysed by |Carrier et al. 
( |20T0l l.~ 



di Mauro et al. (201 1 ) also modelled a red giant observed 



by Kepler that was found in the hydrogen-shell burning phase. 

The hundreds of red giants described by Hekker et al. (2009 ) 
were observed in the so-called "exo-field" of CoRoT (intended 
for a large number of faint stars) but CoRoT observed HR2582 
(HD 50890) in its seismo-field (used for a few bright stars). 
Observations lasted approximately 7obs = 55 days, during the 
first CoRoT observing run dedicated to science. Compared to 
the hundreds of red giants observed in the CoRoT exo-field, 
HR2582 (HD 50890) is an interesting case as its mass can also 
be inferred from non-seismic quantities. It therefore offers the 
opportunity of verifying the compatibility with seismic scaling 
relations. 

HR 2582 (HD 50890, HIP 33243) is listed in the Simbad cat- 
alogue at CDSj^Jas a G6.0III star with an apparent magnitude 
my = 6.03. Some of its characteristics were derived from com- 
plementary observations detailed in Sect. 2.2 which yielded: 
T eff = 4665 + 200 K, [Fe/H] = -0.18 + 0.14, log L/L Q = 2.70 ± 
0.15, and a radius of R/R Q ~ 34 ± 8. Based on the asteroseismic 
scaling relations (Eq.[T]|2|, it appears that HR2582 (HD 50890) 



is more massive than the stars belonging to the red clump men- 
tioned above. This makes this star a rare case since its high 
mass implies a rapid evolution. Modelling this star brings 
insights on the processes at work in the last stages of evolu- 
tion when the star could be burning its hydrogen in a shell 
or its central helium or even the later stage of shell-He burn- 
ing. The brightness and parallax of HR2582 (HD 50890) af- 
ford for relatively tight constraints on its position in the 
Hertzsprung-Russell (H-R) diagram (see Sect. |2.2| >. This is 
not enough to determine if the star is ascending or descend- 
ing the red giant branch (the determination of the observed 
evolutionary stage will require detailed seismic information, 
see Sect.|5) but this nevertheless imposes in turn strong con- 
straints on the modelling of the star's internal structure. This 
is also a good opportunity to assess the uncertainties on the 
results of such modelling. 

Our purpose here is to determine more precisely the charac- 
teristics of the star. We first use the mean large separation (Av) 
and the frequency of maximum oscillation amplitude (v max ) to 
determine the mass and radius of the star from scaling relations 
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Fig. 2. Oversampled (by a factor of 4) power density spectrum 
of the light curve shown in Fig. [T] and a smoothed version in 
thick solid line (using a 1.8-//Hz boxcar). 



together with associated error bars. We study their compatibility 
within the error bars with masses and radii derived from the loca- 
tion of the star in the H-R diagram. These last values also suffer 
from uncertainties in the location in the H-R diagram. In addi- 
tion, they also depend on the physical description used. Thus, 
uncertainties in the physics must be accounted for in the uncer- 
tainties derived for the (non-seismic) mass and radius of the star. 
Despite all of these uncertainties, we show that combining both 
seismic and non-seismic information results in a much more pre- 
cise determination of mass and radius than considering only one 
of these two types of information. In addition, we are able to 
provide a tighter constraint on the metallicity of the star. We also 
derive constraints on the values of the free parameters entering 
the description of convection adopted here. 

Details of the observations and data analysis are given in 
The seismic spectrum is discussed in Sect.[3] A first ap- 



Sect.2.1 



proach to the modelling of the star is discussed in Sect. [4] where 
the influence of the physical description and parameters used in 
the models are estimated, before an optimal model is searched 
for in Sect.|4~2l Age ambiguity for HR2582 (HD 50890) is dis- 



cussed in Sect. 5] Conclusions are given in Sect.|6] 

2. Observations 

2.1. CoRoT photometric data 

A detailed description of the CoRoT instrument and its perfor- 
mance is given by |Auvergne et al.] ( |2009 1 



HR2582 (HD 50890) has been observed for 54.65 days dur- 
ing the first observing run dedicated to science at the beginning 
of 2007. We use in this paper the data ready for scientific anal- 
ysi^J Outliers and other missing data points were replaced by 
interpolated points. The influence of interpolation was checked 
to be negligible. The duty cycle before interpolation was 90%. 



1 Centre de Donnees de Strasbourg, 
www. simbad. u-strasbg . fr/simbad/ 



2 The "HELREG level 2" data corresponding to light curves that have 
been corrected for known instrumental effects and resampled onto a 
regular cadence in the heliocentric frame with one measurement every 
32 seconds ( Samadi et al. 2007 ) 
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Fig. 3. Autocorrelation of the power density spectrum between 9 
and 24 //Hz in solid line (and using an oversampled spectrum 
in dashed line) 



The light curve is shown in Fig. [T] We notice in Fig. [T]a clear 
inverse V-shape that is not explained by any known instrumen- 
tal effect and that is not seen in the light curves obtained for 
other stars during the same observing run. This inverse V-shape, 
whatever is its origin, is a long-term trend that can be removed 
without consequence to our analysis since we are interested in 
periods of less than 1 day. Thus, the data are high-pass filtered 
by substracting two different straight lines, one for the left part 
and one for the right part of the V-shape. 

The power density spectrum of the time series, computed us- 
ing a Fast Fourier Transform, is shown in Fig. [2] After oversam- 
pling by a factor of 4, it has a frequency resolution of 0.05/iHz 
and it is normalised such that the total power integrated from 
zero to twice the Nyquist frequency is equal to the variance of 
the residual light curve; or in other words, the Fourier transform 
is normalised by 1 / Vr where T is the total length of the obser- 
vations. 



2.2. Complementary observations 

The Hipparcos parallax of HR2582 (HD 50890) is n = 2.99 + 
0.44 mas according to the new determination of van Leeuwen 

To determine the atmospheric parameters of HR2582 



(HD 50890), we analysed a spectrum from CORALIE ( Queloz 
|et al.| 2000) obtained on 20 November 2001, which is available 
through the GAUDI database (Solano et al. 2005). The spectrum 
has a resolution of 45000 and signal-to-noise ratio in the contin- 
uum around 6000 A of 160. The lines are relatively broad and 
we measure v sin i = 10 ± 2 km/s. The broad lines made it diffi- 
cult to do the analysis due to blending, but our results are based 
on the calculation of synthetic spectra that take the blending into 
account: we used the VWA software (Brunt fet al.|2010b|a| ). 

We used 44 Fe i (neutral iron) lines while only 2 Fe n lines 
were available. To constrain logg one often uses the condition 
that Fei and Fen abundances must agree. However, this was not 
possible since so few lines of Fe n were available. However, this 
approach was not very robust since only three partially blended 
Fen were available (6084.1, 6247.6, 65 17.0 A). Also, the wide 
Ca lines at 6122 and 6162 A could not be used to further con- 



Fig. 4. Power spectrum of the light curve shown in Fig. [I] once 
the contribution of the background noise has been removed (see 
text). No oversampling is used in order to compute the statis- 
tical levels of confidence. The dashed lines indicate the 99% and 
90% levels of confidence for detection (see text). The crosses in- 
dicate frequencies that are equally spaced by the value found 
from the power autocorrelation, 1.7 //Hz. 



strain log g, due to their lack of sensitivity. With this lack of con- 
straints, we found r eff = 4665 ± 200, logg = 1.4 ± 0.3 and 
[Fe/H] = -0.18 ±0.14. 

On the basis of these data, we estimated a bolometric correc- 



tion of -0.40±0.12 mag following [VandenBerg & Clem| (2003 ). 
It yields a luminosity of log L/L Q = 2.70 + 0.15 and, through the 
Stefan-Boltzmann law, a radius of R/R Q ~ 34 + 8. 

3. Interpretation of the seismic spectrum 

In the power density spectrum, an excess of power is clearly seen 
around v max = (15 + 1)//Hz, showing a bell-shape, enhanced in 
the smoothed version of the power spectrum, a well-known prop- 
erty of a solar-like oscillation spectrum. The central frequency 
v max of this observed excess of power can be used to determine a 
first estimate of the mass of the star, based on the scaling relation 
suggested by |Brown ( 1991| l and derived by Kjeldsen & Bedding 
(1995): 



M/M G 



(R/Re) 2 VTW/Tke 



(1) 



From the stellar parameters of HR 2582 (HD 50890) given in 
Sect.[T] we derive a first estimate of the mass of approximately 
5.2 + 2.9M Q . 

We then derive the statistical significance of the peaks observed 
in the power spectrum (without oversampling, thus giving 
a resolution of 0.21//Hz), using the so-called null hypothesis 
(Fish er|1925 1. The spectrum is considered to be made of white 
noise and a confidence level is computed from the probability of 
a peak caused by noise to reach a given level of power. To apply 
this, the contribution of the background noise in the power den- 
sity spectrum between 9 and 24 //Hz was removed by fitting a 
power law of the frequency and dividing the observed raw spec- 
trum by this. The resulting power density spectrum is shown in 
Fig. [4] The noise level was determined from the mean power in 
the frequency range considered. Then, the significance threshold 
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Table 1. Frequency of the p modes identified using their regular 
frequency spacing. 



Frequency 
QiHz) 
9.8 
11.4 
12.9 
14.6 
16.4 
18.2 
19.8 



Table 2. The uncertainty in frequency is conservatively consid- 
ered to be the frequency resolution (6v = 0.2//Hz). 
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Fig. 5. Echelle diagram of the observed modes (in red) with Av = 
1.70//Hz and the modes computed from the optimal model (in 
black) with Av = 1.72/iHz (see Sect.|4~5]l. 

corresponding to a probability of 90 and 99% for the peaks not to 
be caused by white noise was computed and are represented also 
in Fig. [4] It appears that the highest peak of the spectrum, around 
15 //Hz, is the only one above the 99% level of confidence. 

Another expected characteristic of a solar-like oscillation 
spectrum is its comb-like structure. Thus, we searched for reg- 
ularly spaced frequencies in the power spectrum by computing 
the autocorrelation of this spectrum (Fig. [3]). This autocorrela- 
tion clearly shows an unambiguous signature of regularly spaced 
peaks with a mean spacing value of (Av) — 1.7//Hz. This ob- 
served value can be compared with the expected value derived 
from the well known scaling relation and the first mass estimate 
from Eq.[T] 



<Av) 
<Av > 



M\ 1,2 (R'- 3 ' 2 



, Mr 



<Av> ~ (1.7 ±0.4) //Hz (2) 



which is in very good agreement with what is observed. 

Then, coming back to the power spectrum, one can see that 
regularly spaced peaks follow the 1.7 //Hz spacing (indicated 
in Fig. [4]). However, they have a low S/N ratio, between 1.2 
and 2.3. Thus, from the simple null hypothesis, they cannot be 
considered as being modes, but the null hypothesis is designed 
to evaluate the confidence of detection of peaks only based on 
their S/N ratio. In the present case, the regular pattern that these 
peaks show is a strong argument to consider them as possible p 
modes. However, we considered that their frequencies were not 
very reliable because of this very low S/N ratio. Thus, the fre- 
quencies were extracted simply from the peak maxima in the 
power spectrum and are given in Table [2] and we considered 
the spectrum resolution (5v = 0.2//Hz) to be the uncertainty of 
the frequencies. The individual frequencies will not be used for 
modelling of the star, which will use the mean large frequency 
separation instead. From the frequency list, we can recompute 
the mean large frequency separation and its standard deviation: 
(Av) = (1.7±0.1)//Hz. 

This result is obtained assuming that the observed peaks cor- 
respond to modes with the same degree £, For reasons of mode 
visibility, the most probable degrees to be detected are £ — and 



£ = 1. It seems improbable that £ — are not visible if £ — 1 
modes are. However, in the present case, assuming the presence 
of ( = 1 modes would give a frequency separation twice as 
large and thus not agreeing with Eq.|2] Hence, we deduce that 
non-radial modes are not visible in the present observations of 
HR25 82 (HD 50890). |DeRidder et al.|pX)9) and |Bedding etaL 



(2010 1 have shown that red giants generally show non-radial 
modes in seismic spectra similar to that of the Sun (as also ex- 
pected by Dup ret et al.||2009] for stars of 2 M Q ). In the present 
case, the relative shortness of the observations (55 days) may be 
at the origin of a too low resolution to detect non-radial modes 
(in particular mixed £=\ modes), or the relatively high mass of 
HR2582 (HD 50890) leads to small amplitudes for non-radial 
modes. The low visibility of £ — 1 modes for some red giants is 
confirmed by Kepler observations (see Moss eTet al.|2011c] ). 

In addition, as observed by |Hekker et al.| ( |2009[ ) and dis- 
cussed by Stel lo et al.| ( |2009| l, the frequency spacing Av and the 



frequency of the maximum of the spectrum, v max , are linked. 



The comparison with Fig. 6 of |Hekker et al.| ( |2009 l shows that 
the observed value of v max is in agreement with the frequency 
spacing observed if one considers that only radial modes are vis- 
ible. More precisely, the Av observed for HR2582 (HD 50890) 
corresponds to the lower values computed from the correlation 
observed by Hekker et al. (2009). This is another indication of 
the high mass of HR2582 (HD 50890) be cause the relation be - 
tween v m .„ and Av is sensitive to the mass (Hek ker et al.|201 1) . 



By comparing 



and Av, one can estimate that this star 



oscillates with low-order n modes (approximately 5 < n < 12), 



which will be confirmed by modelling in Sect. 4.5 



We then used the mean value of the large spacing to cut the 
power spectrum into pieces of 1.7/zHz. These pieces are then 
stacked on top of each other to build an echelle diagram, shown 
in Fig. [5] where modes with same degree £ appear in a ridge more 
or less vertically aligned. However, this ridge has a typical cur- 
vature that is observed in other stars, in particular in the Sun or 
solar-like pulsators and that will be compared with expectations 
drawn from stellar modelling in the following section. 

4. Stellar characteristics for HR2582 (HD 50890) 

The modelling of the star is carried out using CESAM2k ( |Morel| 



1997 Morel & Lebreton 2008 1. In order to compute a stel 



lar model for a given mass and age (or luminosity and T e ff), 
one needs to specify its initial relative helium abundance Y 
and its metallicity Z/X. One must also provide values for the 
free parameters used in the physical description of the stellar 
interior. In the present work, convection is included using the 
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Fig. 6. Left: Evolutionary tracks for models with masses from 2.8 to 5.5 M Q and observational luminosity and effective temperature 
constraints (blue box); Right: zoom of left panel with models satisfying the seismic constraints shown with coloured points (cyan 
indicates models having a compatible v max , magenta a compatible Av, and dark green models satisfying both conditions). 



Canuto, Goldman, Mazzitelli (CGM) formulation (Canut o et al 
1996). This formulation requires a prescription for the charac- 
teristic scale length associated with the energy-bearing eddies. 
Following Bernkopf|( |1998| l, we assume that this scale length is 
Q'cgm H p where H p is the pressure scale height and q'cgm a free 
parameter. 

Another free parameter is needed to describe convective 
overshoot since stars as massive as HR2582 (HD 50890) have 
a convective core while on the main sequence. The modifica- 
tion of the luminosity when convective core overshoot is in- 
cluded during the main sequence is maintained once core con- 
vection has disappeared after the main sequence. Hence core 
overshoot must be taken into account in modelling the red gi- 
ant HR2582 (HD 50890). Overshooting is described as an ex- 
tension of r c , the radius of the core as defined according to the 
instability criterion of Schwarzschild, over a distance taken to be 
d ov - a m mm(r c ,H p ) where a ov , the overshooting parameter, is 
our second free parameter. 

The models considered are computed with the OPAL 2005 
equation of state and opacity tables as described in |Lebreton 



(see Sect. 4.1 1. Then, the influence of the chemical composition 



et al.| ( 2008| l. The chemical mixture by [Asplund| ( 2005) > is" 



adopted for the calculation of the opacities. For the calculation 
of the equation of state, we have assumed the chemical mixture 
of |Grevesse & Sauval| ( |1998| l with the except ion of CNO for 
which values close to the recommendation of Asplund (2005) 
are adopted. Neither microscopic diffusion nor near-surface ef- 
fects are included in the modelling. 

First, the star is modeled for a given physical description 
(fixed values for a ov , Q'cgm) and a given chemical composition. 
The "best" model is searched in a pre-existing grid of models 
for this description, defined as the model that minimises the dis- 
crepancies with observational constraints, given this description 



and physics parameters on the mass and radius of the best model 
are investigated individually (Sect. 4.2 1. Finally in Sect. 4.2 the 
minimisation of discrepancies between models and observation 
is performed by adjusting several of the free parameters. 



4. 1 . Modelling for a given physical description 

The models considered here are built assuming the same chem- 
ical composition as the Sun, that is the metallicity Z/X = 0.018 
(Z = 0.0135) and helium abundance Y - 0.2492 and the same 
mixture. We also assume q'cgm = 0.763, which is the value 
obtained for the Sun with the physics used here. As HR2582 
(HD 50890) is expected to have a convective core, we further as- 
sume convective core overshoot with o- ov = 0. 1 on the main se- 
quence. Fig. [6] shows that stars with different masses (from ~ 2 
to ~ 5.5M Q ) have evolutionary tracks compatible with the ob- 
servational contraints L and T e s. In this range, higher masses 
correspond to less evolved stars. All of the models describe stars 
on the ascending red-giant branch. However, models describing 
more evolved stars, at the helium-burning stage, can have the 
same position in the H-R diagram. This point will be discussed 
in Sect.[5] 

Seismic information, such as the mean large separation, al- 
lows a more precise characterisation. The mean large separa- 
tion links the mass to the radius through Eq.|2] and thus corre- 
sponds to a given strip in the H-R diagram (see for example 
Fig. [8}. This allows a strong limitation of the possible observa- 
tional values, which is quite useful as the luminosity is poorly 
estimated because of the relatively large uncertainty on the par- 
allax of HR 2582 (HD 50890). We compute the mean large sep- 
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Fig. 7. Computed values of x 2 (Eq[3j for models corresponding to a given physical description as a function of the mass. On the 
left, the constraints are (Av) > and v max , in the middle panel r e ff is added, and on the right the luminosity is also taken into account. 
Note the different scale in the latter case. 



Table 3. Characteristics of the models presenting the lowest x 2 
value in Eq.[3]for a given physics and chemical composition and 
when using a varying number of constraints (see text). 



Constraints 


(Ay) 


(Av), y max 


(Av), v max 


used 




T e ff 


Ten, log L/L, 


2 

X 


0.008 


0.36 


< 1.1 


M/M 


3.6 


3.3 


3.8-4.0 


R/Re 


28.3 


27.5 


29.4 - 30.9 



tainty on the parallax of the star. However, it contributes addi- 
tional information on the radius and thus should be taken into 
account. The consequence is a larger uncertainty on the mass 
determination as the lower values of x 2 can be reached for mod- 
els of different masses. 



4.2. Influence of the physical description on the modelling 



aration (Av) and v max according to Eq.[T] and [2] for our model 
grids. The acceptable models (that is those that satisfy the ob- 
servational constraint on (Av)) correspond to the highest tem- 
peratures and lowest luminosities of the observational range for 
HR2582 (HD 50890) and lead to a better mass estimate: ~ 3 to 
~4M Q . This is confirmed by the additional constraint provided 
by the frequency of the maximum power, v max . This might not 
come as a surprise as these two observables have been shown 
to be correlated for cool giants (Hekk eret al.|200 9; Stello et al. 
[2009] l. 

We now wish to determine which model provides the "best" 
match with all of the observations. For all of the grid models, we 
compute the quantity^- 2 : 



X 



Z 



(3) 



where Xj are observables (obs and mod specifying the ob- 
served or the modelled value) and cr, the associated uncertainty. 
Computations of x 2 were first performed with X\ = (Av) and 
X2 = v max , then adding X3 = T e ff and finally also taking into ac- 
count^ = log L/Lq. The theoretical seismic quantities (Av) and 
v max are computed using Eq.[T] and [2] with the mass, radius and 
effective temperature obtained from the evolutionary calculation. 
The resulting values for^f 2 are listed in Tableland can be seen 
in Fig. [7] The mass estimate is M = (3.9 ± 0.1)M o and the other 
resulting model parameters are listed in Table|3] A given param- 
eter can be both a constraint with an uncertainty (thus having to 
be in a given range), and an output. 

When the luminosity is taken into account, the different x 2 
values are significantly higher and form a broader shape in Figj7] 
low values of x 2 can be reached for several models. The photo- 
metric luminosity constraint is less robust because of the uncer- 



A change in the assumed physical description is likely to modify 
the stellar characterics of the "best" models determined in the 
previous section. We illustrate this issue by assessing the effect 
of changing the values of the free parameters related to convec- 
tive motions, qcgm and a ov , as well as the metallicity Z. 

An illustration of the influence of varying qcgm (with all 
other parameters kept constant) on evolution tracks in the H-R 
diagram is shown in Fig. [8] The models satisfying the obser- 
vational constraints in the H-R diagrams shift to lower lumi- 
nosity and temperature with q'cgm increasing. In order to take 
into account the influence of other parameters, model grids with 
or ov = and or ov = 0.2 were computed, everything else being 
kept the same as for q ov — 0.1, as in Sect. |4.1| Using the same 
approach as above, x 2 (for all four constraints) are computed and 
displayed in Fig.|9]where the effect of the variation of q'cgm on 
mass determination appears clearly: a higher q'cgm leading to a 
lower mass, the variation being as high as 1 .0 M for a variation 
of 0.2 in Q'cgm- The smallest value of qcgm=0.6 provides the 
lowest values of x 2 - Proceeding in the same manner for varia- 
tions of q ov , mass variations are smaller but still of the order of 
0.2 M Q . The different values of q ov lead to similar values of^ 2 . 
The metallicity also has an influence on the mass determination: 
varying the solar metallicity to [Fe/H] = -0.18 the difference 
in mass reaches approximately 0.2 M G . The lowest^ 2 values are 
clearly in favour of the solar metallicity, but even lower values 
could be obtained for intermediate values of the metallicity. 

The different stellar characteristics, yielded by estimating the 
lowest x 2 values in a grid of models for the different parameters 
considered here, are summarized in Table|4] These values do not 
aim at being precise (the precise outputs of the optimum model 
are given in Sect. 4.3 1 but at showing the approximate range of 
the model parameters. 
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Fig. 8. Evolutionary tracks for different masses and different qcgm (left: 0.6; center: 0.7; right: 0.8). H-R constraints (L and T e g) 
are shown with the blue box, and models satisfying the seismic constraints are shown with coloured points (cyan indicates models 
having a compatible v max , magenta a compatible (Av), and orange models satisfying both conditions). Models with higher values of 
Q'cgm shift compatible masses to lower values (see text). 



Table 4. Sensitivity of the stellar characteristics to the model 
parameters. The estimates presented correspond to the models 
giving the lowest^ 2 . 








0.1 


0.2 


M/M Q 


3.9 - 4.2 


3.6-4.0 


3.6-4.0 


R/Re 


29.5 - 30.5 


28-30 


27.5 - 30 



Q'CGM 


0.6 


0.7 


0. 


i 


M/M Q 


4.4-4.8 


3.6-4.1 


3.3- 


4.0 


R/Re 


31-32 


29-30 


27- 


29 



[Fe/H] 





-0.18 


M/M Q 


3.6-4.0 


3.5-3.7 


R/Re 


27.5-31 


27-28 



the mass and from similar results not shown for the radius. They 
do not aim at a quantitative comparison but at an overall view 
only. 



4.3. Optimal modelling 

As shown in the previous section, the search of an optimal model 
should allow all of the parameters related to the physical de- 
scription to vary. However, this can suffer from limitations men- 
tioned below. The determination of the best model used here is 



perfor med with the method described in Miglio & Montalban 
( |2005| l: x 1 i s determined from a minimisation process using a 
Levenberg-Marquardt algorithm coupled to the CESAM algo- 
rithm in order to converge to the parameter values that provide 
the lowest^- 2 . The searched models are restricted to the ascend- 
ing branch of red giants and do not consider the further evolution 
of the star (He-core burning, this will be discussed in Sect. [5}. 

However, the application of this method cannot avoid some 
limitations. First, many parameters are needed to fully describe 
a model whereas there are in comparison only a few independent 
observational constraints. If all parameters are free, the result of 
the minimisation is not reliable as it depends on the initial values 
of the parameters. This led us to run the minimisation process for 
a grid of fixed values of a given parameter (see below). Second, 
there is some cross-talk between the model parameters (for ex- 



ample between the mass and Y, and also between Q'cgm and Y). 
The consequence is a degeneracy of the solution. 

A first set of minimisation processes was run with different 
values of [Fe/H]. It appears that extreme values of the metal- 
licity within the observational range lead to values of Y either 
lower than the primordial He abundance ([Fe/H] = -0.18 yields 
Y = 0.24), or too high ([Fe/H] = yields Y = 0.30, a value that 
corresponds to those of the youngest open cluster stars and there- 
fore lies in the upper part of the reasonable range). Consequently, 
the minimisation will be then processed using a fixed, interme- 
diate value of the metallicity: [Fe/H] = —0.1. 
Then, the minimisation process was run with fixed values of the 
metallicity, Y, Q'cgm and q ov which could take three values: 0, 
0.1 and 0.2. The best values of^f 2 were obtained for a m =0.1. 
The next parameter to be investigated was Q'cgm- A series of 
minimisation processes was run for decreasing values of qcgm 
(from 0.75 to 0.44, but fixed for each minimisation process run) 
whereas M, Y and the core temperature (from which the age is 
derived) are free parameters. A clear trend appears when select- 
ing the models giving the lowest values of^f 2 (smaller than 0.3): 
they correspond to low values of q'cgm (smaller than 0.6, see 
Fig. 10 1. In addition, when selecting these models (y 2 < 0.3), a 
clear anti-correlation appears between their M and Y values (see 
Fig.[TT}. 

A similar relation is observed between q'cgm and Z. For the opti- 
mum model, Qcgm — 0.45 ±0.15. It is obtained for a fixed metal- 
licity: [Fe/H] = -0.1. If the metallicity is fixed at [Fe/H] = 0, 
then qcgm - 0.51 +0.15 (see fig. 



10) 



Both values are substan- 
tially lower than the solar value (qcgm.q — 0.76) but are fully 
compatible with results from 3D numerical simulations of con- 
vection in red giants (qcgm.g - 0.6, Samadi et al. - in prep.) 

The characteristics of the optimal model (lowest x 2 ) are 
listed in Table|6j but one should not forget that this model corre- 
sponds to given values of metallicity and q ov . 

(Av) has been computed from the mode frequencies given by 
the optimum model, and an expected v max is then computed us- 
ing the scaling relation between these two parameters. As seen 
from Table[6] the computed v max is in agreement with the ob- 
served value. 

A last remark on this optimal model: it was found with a 
metallicity in the upper part of the interval given by observa- 
tions. As mentioned in Sect.|2.2| the parameters T e s and [Fe/H] 
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Fig. 9. Computed values of x 1 (Eq.[3]> versus the mass, when one 
parameter of the modelling is varied while the others are kept 
constant: acGM (top, M/M Q = 4.6 ± 0.2), or ov (middle, M/M Q = 
3.9 ± 0.3) and [Fe/H] (bottom, M/M Q = 3.8 ± 0.2). The final 
estimate is 3.6M Q < M < 4.8M . 



Fig. 10. Variation of x 1 in the minimisation process for differ- 
ent fixed values of acGM (and a ov =0.1, [Fe/H] = (in blue) 
and -0. 1 (in red) fixed, whereas the mass, age and initial helium 
abundance are free parameters). 



Table 6. Characteristics of the optimal model (giving^- 2 = 9.35 ■ 
10~ 2 ). Uncertainties are listed in Table^ 



Constraints 


log(L/L ) 


2.70 ±0.15 


iog(r cff ) 


3.67 ± 0.02 


Vmax 0"HZ) 


15.0 ± 1.0 


Av (jiRz) 


1.7 ± 0.1 


Fixed parameters 


[Fe/H] 


-0.1 




0.1 


a CGM 


0.46 


Fitted parameters 


M/M Q 


4.63 ± 0.09 


Y 


0.298 ± 0.034 


Tcore (MK) 


80.3 ± 14 


Outputs 


of the model 


Age (Myr) 


105.5 


log(L/L ) 


2.66 


iog(r cff ) 


3.67 




15.0 


Av (pHz) 


1.69 


R/Re 


32.27 



were determined without strong constraints on logg. Using the 
values of M and R from Table|4] it was then possible to itera- 
tively improve the atmospheric parameters using log g from the 
seismic modelling, which is close to 2.09 ± 0.02. For a fixed 
value of logg, we found r eff = 4700 ± 180K and [Fe/H] = 
-0.05 ±0.11. This higher observed metallicity obtained a poste- 
riori is in agreement with seismic modelling. 
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Fig. 11. Values of Y versus M for models obtained with the min- 
imisation processes (with fixed values of (Xcgm) and closest to 
the constraints (x 2 < 0.3): a clear anti-correlation appears. 

4.4. Discussion of the accuracy of HR 2582 (HD 50890) 
stellar parameter determination 

As discussed in the previous sections, there are several ways 
to estimate the mass of a star, and this determination can rely 
on different observational constraints. We have first used model 
grids. The corresponding evolution tracks in the H-R diagram 
can be compared to the observed luminosity and temperature 
and their error bars. Then, seismic constraints (v max and Av) can 
also be taken into account to select the adequate evolution tracks 
in the H-R diagram. This selection can be quantified by look- 
ing at the sum of squared differences x 2 between observed and 



computed parameters from a grid of models (see Sect. 4.1 1. This 
could be performed with a grid of models computed for a given 
"physics description" (parameters describing the overshoot, the 
convection and also a given metallicity and helium abundance) 



or allowing this physics to vary (see Sect. 4.2 1. Finally, an op 



timal modelling can be performed, allowing the physics param- 
eters to vary during the process of the minimisation of x 2 (see 
Sect. 4.3 1. To be complete, we recall that the large separation as 
an output from modelling can result from the use of the scal- 
ing relation (Eq.[2| and the modeled effective temperature, or 
directly from the computed modeled frequencies. 

We could have also used the seismic constraints only, to- 
gether with the observed temperature and combining the scaling 
relations of Eq.[T]and[2] Then, one obtains the following uncer- 
tainties on the mass determination: 



QjM) = g gtW) | ^(Av) 3 ojTes) 
M v max Av 2 r eff 

and similarly for the radius: 

cr(R) cr(v max ) | 2 cr(Av) 1 o-(T eB ) 



R 



Av 



2 r, 



(4) 



(5) 



eff 



However, these scaling relations are empirical and hence suffer 
from some uncertainty. The corresponding estimates are given 



Fig. 12. Upper panel, solid line: Av sc , computed from the scaling 
relation (Eq.[2]i and the modeled mass and r e jf versus Av com- 
puted from the actual frequencies given by the optimal model 
(the blue frame indicates the observations); lower panel: Av sc 
and Av versus the computed v max from Eq.[T]for a series of mod- 
els satisfying the observations represented as the blue frame. 
This validates the use of the scaling relations in the computa- 
tions. 



in Table|7] together with all of the other resulting mass esti- 
mates and their uncertainties. One can draw some conclusions 
and some limitations from this table. As expected, the use of the 
seismic constraints together with the temperature and luminosity 
substantially improves the mass determination. However, seis- 
mic constraints alone (using the scaling relations) do not yield 
a precise estimate. Such a precise estimated is given by the use 
of a grid of models, which can however lead to a largely under- 
estimated uncertainty if the model physics is fixed. 



4.5. Echelle diagram 

The echelle-diagram built with the frequencies of the optimal 
model is compared to the observed one in Fig. [5] The theoretical 
echelle-diagram (computed without corrections of near-surface 
effects) matches the observed one well within the observational 
error bars. Note that the folding value for the model, 1.715 //Hz, 
is only slightly different from that of the observations (1.7yuHz). 
This value, usually associated with the mean large separation 
is different from the value 1.7/vHz given by Eq. [2] which is the 
relation used in the minimisation process. It nevertheless falls 
within the uncertainties associated with this relation. The present 
modelling also confirms the observed values of radial order (5 < 
n < 12). 

We also checked that the estimation of Av from a scaling 
relation and the modeled mass (Eq.|2]i was correct by comparing 
with the frequencies derived from the model (see Fig.[T2|. 
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Table 7. Summary of mass estimation and uncertainty. Av represents the large separation derived from Eq|2| (Av) represents the 
large separation computed from the modeled frequencies. T°^ s and T™ d are respectively the observed and modeled temperatures. 



Method 



Parameters used 



M cr(M) cr(M)/M(%) Mass interval (M ) 



R cr(R) cr(R)/R(%) Radius interval (R e ) 



Model grids and 
H-R diagram 



yobs yobs 



3.5 
34 



57 
24 



1.5 < M < 5.5 
26 < R < 42 



Scaling 



(using Eq 



and 



with r£ s ) 



Av, v max 



3.6 

28.8 



1.8 
6 



50 
21 



1.8<M<5 
22.6 < R < 34.6 



(using Eq 



]\£qdel sad with fixed input parameters 
2 



and 



with r e 7 d , 



TTnod 
J eff 



compatible with T°g") 



T° hs Av v 



3.9 

29.2 



0.1 
0.3 



2.5 
1.0 



3.8 < M < 4.0 
29.0 < R < 29.5 



Medel grid with varying input parameters 



(using Eq 



and 



with T™, r c 7 d compatible with r°° s ) 



robs •robs A 

' eff ' ' m 



4.2 
30.0 



0.6 
1.4 



14 
4.7 



3.6 < M < 4.8 
28.5 <R < 31.3 



-i Optimal modelling 



(using Eq 



with T™\ r™ d compatible with r* s ) 



robs jobs A 

' eff ' ' max 



4.2 



0.9 



29.9 1.9 



21.5 
6.3 



3.3 < M < 5.3 
28.0 <R< 31.8 



5. Evolutionary stage ambiguity 



As mentioned earlier, the location of HR 2582 (HD 50890) in 
the H-R diagram can correspond to several evolutionary stages: 
either on the H-shell burning on the first ascending branch or He- 
core burning on the descending or second ascending branch (see 
Fig. 13 i. The respective ages for these three evolutionary stage 
are 157, 163 and 180 Myr. The corresponding helium abundance 
in the core are Y c = 0.983 (no helium fusion yet), 0.768 and 
(helium fusion phase ended). At a given location the star has 
similar mass and radius and the large separation Av is unable 
to discriminate between the possible evolutionary stages, as dis- 
cussed by Mazumdar et al. ( 2009| l. To emphasise this problem 
in the present case, we have considered a 4 M Q star computed 
assuming a CC M = 0.7, a m = 0.15 [Fe/H] = 0, Z = 0-0172 , 
Y = 0.28 using the solar mixture of |Grevesse & Noels| ( |1993} . 
The values of (Av) (computed as the mean separation of modes 
with 5 < n < 12) for the three stages of evolution are 1.53, 
1.49 and 1.52yuHz, falling in a quite narrow range, smaller than 
the frequency resolution of the present data set (6v = 0.21yuHz). 
The three models do not have exactly the same radius, hence not 
exactly the same mean large separation. This explains why the 
I — ridges in the echelle diagram shown in Fig.[l4]are slightly 
shifted from each other. If one scales the offset with the respec- 
tive mean large separation for the three cases, the variations of 
the three ridges with the frequency nearly coincide. This con- 
firms the expectation that the variations of the ridges with the 
frequency arise from the properties of the surface layer which 
are similar for the three models (see Fig.[T4|>. 
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3.7 3.693.68 



4 3.8 
log T ef[ 



Fig. 13. Evolutionary track for a 4 M Q model showing that it re- 
turns to the same location of the H-R diagram at the beginning of 
the ascending branch (lowest track) and during He-core burning 
(descending - intermediate track - and second ascending (higher 
track) branches). Selected models discussed in the text are plot- 
ted with solid black dots. 



If non-radial mixed modes were present, and if the time 
series were long enough to provide a sufficient frequency res- 
olution to resolve period spacings of these mixed modes, the 
age ambiguity could be removed thanks to the age signature in 
mixed-mode spacings as shown by Beddi ng" et al.| ( |201 l| l and 
Moss er et al.| ( |20lTa] ). However, for a star with such a low (Av) 
value, this requires an observing time longer than 2 years. 



6. Conclusions 

HR2582 (HD 50890) is a luminous giant that shows a number 
of peaks in the power spectrum at low frequency, centred around 
the frequency v max = 15 /^Hz, with a mean frequency separa- 
tion (Av) =1.7 //Hz. This separation can be expected from scal- 
ing relations (Eq.[T| and |2]i and implies that only radial modes 
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Fig. 14. Echelle diagram for the 4 M models of Fig 13 (Av) is 
taken to be 1.53, 1.49 and 1.52 //Hz (for the observed range: 5 < 
n < 12) for the three stages of evolution: the youngest (violet), 
intermediate (red) and the oldest (black) ones. The differences 
between the three cases are too small to be detected with the 
present data set. 



are cle arly observed. De Ridder et al. (2009 1 and Bedding et al. 
( 2010| > clearly showed the presence of non-radial modes in red 
giants. Their absence in our observations could be caused by a 
too short observation sequence or, since our results for HR 2582 
(HD 50890) lead to an estimated mass around 4 M Q , the internal 
structure (namely the ratio of the density of the core to the mean 
density) could be different and lead to different amplitudes for 
non-radial modes in high-mass stars. 

The modelling of the star is shown to be sensitive to the phys- 
ical description used (parameters describing convection, over- 
shoot, etc..) and other constraints such as the metallicity for ex- 
ample. For example, the value of the estimated mass can vary by 
about 1 M for different parametrization of the convection. Only 
some aspects of the physics involved in stellar modelling were 
varied in the present work. Thus, the uncertainties derived from 
modelling are possibly under-estimated: for example, the influ- 
ence of microscopic diffusion was not taken into account here. 

An optimal model was found with a mass of (4.15±0.85)M Q . 
Supposing the star is on the ascending red giant branch, its age 
was estimated from this optimal model to be 155 Myr. Another 
interesting result is the low value of o-cgm. lower than both the 
solar value and results obtained from simulations, indicating 
that the modelling of red giants could bring in return particular 
inputs on the description of convection. 
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The echelle-diagram built from the optimal model is globally 
in agreement with the observed one (Fig. [5]). The small differ- 
ences between modeled and observed echelle-diagrams contain 
potentially more information about the internal structure of the 
star and will be investigated in a forthcoming work. For exam- 
ple, microscopic diffusion is known to lead to strong element 
depletion for such massive stars. Then, turbulent mixing should 
be present to compensate for such a non-observed high element 
depletion. It is possible that seismology of high mass stars such 
as HR 2582 (HD 50890) can bring some constraints on this issue 
(to be investigated in a future work). 



